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Abstract 

Based on numerical simulations, we show that higher-order nonlinear indices (up to ns and 
nio, respectively) of air and argon have a dominant contribution to both focusing and defocusing 
in the self-guiding of ultrashort laser pulses over most of the spectrum. Plasma generation and 
filamentation are therefore decoupled. As a consequence, ultraviolet wavelength may not be the 
optimal wavelengths for applications requiring to maximize ionization. 
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The filament at ion of high-power, ultrashort laser pulses is generally described as a dy- 
namic balance between Kerr self-focusing and the defocusing by the plasma generated at 
the non-linear focus jl-5]. In this description, the Kerr-induced change in the refractive 
index is truncated to the first term n 2 I, where I is the local intensity. However, at higher 
intensities, the development has to be extended to higher-order terms in J, so that the real 
part of the refractive index at any frequency u writes: 



n{tu\ 
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where the n 2 j(co) coefficients are related to the (2j+l) th order susceptibility tensor x( 2j+1 )(u;), 
in the degenerate case where all considered fields are at frequency u: 
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n ( w ) = .^V - Me (x (2j+1) (^)) (2) 

In the last years, several numerical works have investigated the influence of the quintic 
nonlinear response on the filamentation dynamics at a wavelength of 800 nm, although 
without knowledge of its actual value jfj-lO]. They suggested that this term was defocusing, 
but considered it as marginal. Recently, the measurement of the higher-order Kerr indices 
at a wavelength of 800 nm up to rig in air and riio in argon [ll| showed that they have 
alternate signs, and are therefore alternatively focusing and defocusing. Furthermore, by 
implementing these terms in a numerical simulation of filamentation, we have recently shown 
that the defocusing terms and rig, rather than the plasma, provide the main regularizing 
process in the filamentation of 30 fs pulses in air at 800 nm, so that plasma generation 



and propagation equations are almost decoupled 



to measurements of plasma-free filamentation 13 



2 [.This finding provided an explanation 
14j and predicts symmetrical temporal 



pulse shapes, in contrast with a balance between the instantaneous Kerr term and the time 



integrated plasma contribution, which implies strongly asymmetric pulse shapes [15|. It 
also explained the discrepancy, by almost two orders of magnitude, between experimentally 
measured and numeric ally derived electron densities within filaments. While the experiments 
yield some 10 14 cm -3 
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171 1. the numerical simulations require a few 10 16 cm 3 to balance 
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the n 2 I Kerr self- focusing term |3|, |4J. The consideration of the higher-order Kerr terms also 
turned out to be necessary to obtain a quantitative agreement between numerical simulations 



and experimental results a 



filled hollow-core fiber 
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aout the propagation of ultrashort infrared pulses in an argon- 
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However, the ionization rates are higher at shorter wavelengths, so that ionization is 



generally believed to be much stronger in the case of ultraviolet filamentation 



20|, while its 



low efficiency in the infrared would expectedly prevent self-guiding of filaments. But the 



recently derived generalized Miller formulae 



21] predict larger absolute values of the higher- 



order non-linear indices at shorter wavelengths, with faster spectral dependencies for the 
higher orders. As a consequence, the relative contribution to self-guiding cannot be easily 
extrapolated from qualitative discussions. 

In this Letter, we numerically investigate these relative contributions from the near- 
ultraviolet to the near-infrared. Using the values of the n^j indices at any wavelength as 



obtained from the generalized Miller formulae 2l|, and the values recently measured at 



800 nm 111 ], we simulate the filamentation of laser pulses from 300 nm to 1600 nm in air 
and argon. We show that filaments are efficiently generated at all wavelengths, even in the 
infrared. Moreover, the plasma marginally contributes to self-guiding in the filamentation 
of 30 fs laser pulses in argon for wavelengths longer than typically 400 nm. In air, where 
the ionization of oxygen is about 100 times more efficient than in argon, the Kerr terms 
still provide the largest defocusing contribution, although the plasma contribution can be 
neglected over the whole propagation length only in the infrared. This finding strongly 
impacts the plasma generation. We also find that ultraviolet wavelengths are not optimal 
to maximize ionization in a filament. 



Numerical simulations are performed, as described in detail in 12j, by solving a non- 
linear Schrodinger equation taking the higher-order Kerr terms into account. The reduced 
scalar envelope e defined such that \e(r, z, t)\ 2 = I(r,z,t), I being the intensity, is assumed 
to vary slowly in time and along z: 
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where u and k=rtQu/c are the angular frequency and wavenumber of the carrier wave, k" 
accounts for the linear group-velocity dispersion, p is the electron density, p Qte the density 
of neutral molecules of species £, p c = m e e u!/ql the critical plasma density and m e and 



q e the mass and charge of the electron, no is the linear refractive index at A, W^(|e| 2 ) is 
the photoionization rate of species I with ionization potential Ue, <jg is the cross-section 
for avalanche ionization as defined below in Equation (|S]) and t refers to the retarded time 
in the reference frame of the pulse. The right-hand terms of Eq.Q account for spatial 
diffraction, second order group- velocity dispersion (GVD), Kerr self-focusing (including the 
self-steepening term), defocusing by the higher-order nonlinear refractive indices, plasma 
defocusing, inverse Bremsstrahlung and multiphoton absorption respectively. We neglect 
the delayed orientational response, which for pulses longer than 100 fs would increase the 
self-focusing term j^, 4] and affect the ionization efficiency of N 2 and O2 molecules by less 



than 20 and 10%, respectively 22j, |23j. These opposite effects are negligible in the numerical 



simulations for 30 fs pulse duration, and do not affect qualitatively the discussion below even 
for longer pulses. We also neglect space-time focusing. However, in the following we mainly 
focus on the self-guiding process, peak intensity and ionization, which are little affected by 
this approximation. 

The spectral dependence of the NLSE ([3]) stems from the plasma contribution as well as 
the dispersion of the linear (no) and non-linear (n 2 j) refractive indices. The latter can be 
deduced at any frequency uj from those of O2 and N 2 at the same frequency by following 
the Lorenz-Lorentz relation, which links the refractive index of a mix of non-polar gases to 
its components 241 ]: 

where n(u) is defined in Equation (TjQ) and pi denotes the relative abundance of the species % 
in the mix. Considering that An^err Cl, a first-order development of Equation (T4]) yields, 
for any order j: 

n 2 j,Air(u) = 0.79 • n 2j ,N 2 (uj) + 0.21 • n 2j ,o 2 (w) (5) 
The values of the n 2 j for argon, nitrogen and oxygen can be calculated at any frequency 



u from their values measured at An=800 nm 



and the dispersion curves in thoses gases 



251 ] . through the generalized Miller formula 21] expressed in term of non-linear refractive 
indices: 

, , / n|(w) - 1 \ 2<J+1) 

While n 8jJ 4j r varies by a factor of 2 between 300 and 1600 nm, ionization displays an even 
much steeper spectral dependence. More specifically, the temporal evolution of the electron 
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density is given by 



dt ~ ^ 



=0 2 ,N 2 
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where po,nt £ is the initial density of neutral molecules of species I. Here, attachment to 
neutral molecules and recombination with positive ions have been neglected owing to the 
short pulse durations considered in this work. The cross-section for avalanche ionization is 
calculated on the b aS1S of Dntde, theory S: 
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u 2 m e e n (u)c 

where v e ^ is the mean collisional frequency of an electron with the species i (i.e., the average 
electron velocity, divided by the mean free path of an electron, assuming that only species I 
is present), and c is the speed of light. At atmospheric pressure in air, v Cy o 2 = 1/(1.75 ps), 
^ e ,N 2 = 1/(440 fs), and v Ct A_? = 1/(350 fs). The multiphoton/tunnel ionization rates are given 
by the multi-species generalized Keldysh-Perelomov, Popov, Terent'ev (PPT) formulation 
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where, expressed in atomic units, E p 
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and rri£ are the orbital momentum and the magnetic quantum number, respectively. In air, 
It = rri£ = [26 1 . Zi is the residual ion charge accounting for the difference between the O2 
and N2 molecules and their atomic counterparts. These empirical coefficients Zq 2 = 0.53 
and Zf<s 2 = 0.9 have been measured at 800 nm j^, and are expected to be independent from 
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wavelength. Since argon is an atomic gas, Za t = 1- The factors IC^*^*] and f(l^mi) are 
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To quantify the relative contributions of the higher-order Kerr terms (n^ through nio) 
and of the ionization to the defocusing which balances the self-focusing, we define the instan- 
taneous ratio £ between the refractive index changes induced by both of these contributions 
at any location r: 

im = \Y,n 2j im j \i p ^- (is) 

This expression is in fact the ratio of the magnitude of the two defocusing terms in the 
propagation equation ((3]), namely those accounting for higher-order Kerr terms and plasma 
defocusing, respectively. When considering the overall action of both effects on the whole 
pulse duration, we define a pulse-integrated value of £: 

^ = / E^/^'l • \ £ (r,t)\dt/ J ^ \e{?,t)\dt, (14) 

We numerically integrated the propagation equation (J3J) for an ultrashort pulse typical 
of laboratory-scale filamentation experiments: 30 fs Fourier-limited FWHM pulse duration, 
a peak power of 6.5 critical powers Per = A 2 /47rn2, a beam diameter of or = 4 mm, a focal 
length f = 1 m and a pressure of 1 bar of either air or argon. The main results are displayed 
in Figure [TJ 



Our simulations yield filamentation over the whole investigatec 



spectral range (300 - 1600 



nm), reproducing experimental observations from the UV [28|, |29] to the mid-IR 30] . As is 
visible in Figured] (a) and (b), the clamping intensity, filament onset and filament length are 
very similar in air and argon over the whole considered spectral length for a given incident 
reduced power P/P cr . Moreover, longer wavelengths yield longer filaments with an earlier 
onset. Furthermore, the output spectrum after 1.5 m propagation (panels (e) and (f)) is 
broader for longer wavelengths. As a result, the relative broadening, defined as the ratio 
of the output spectral width to the initial frequency, Au/uq, is almost constant across the 
spectrum. Spectral properties appear very similar in air and argon. 

6 




FIG. 1: Maximum on-axis intensity (a,b), electron density (c,d), spectral broadening (e,f) and 
pulse-averaged ratio £ of the on-axis pulse-averaged contributions of the higher-order Kerr and 
plasma contributions to the non-linear refractive index (g,h; see text for details), as a function of 
the wavelength and propagation distance in air (a,c,e,g) and argon (b,d,f,h). 

In contrast, the electron density is approximately 10 times higher in the former than in 
the latter. The spectral dependence of the electron density is very non-monotonic (Panels 
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FIG. 2: Wavelength dependence close to 800 nm of the electron density generated by a 30 fs pulse 
of constant intensity 50 TW 
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FIG. 3: Wavelength dependence of the clamping intensity canceling the Kerr effect /An Kc rr=0; m 
0*2, N2, air and Ar, based on the dispersion of the higher-order Kerr indices from Equation ([6]) 

(c) and (d)). On one hand, due to steps in the number of photons required for ionization, 
the electron density is not a monotonic function of the wavelength. For instance, ionization 
in air is almost 2 times more efficient at 793 nm than at 815 nm (See Figure [2]). On the other 
hand, contrary to expectations that shorter wavelengths should result in stronger ionization, 
a maximum in the peak electron density is observed in the 610 nm range in the case of air, 
and around 470 nm in argon (Figures [T](c),(d) and|U(a)). Away from this maximum, the 
peak electron density decreases within a dynamics of typically one order of magnitude. This 
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unexpected behavior stems from the convolution of the respective spectral dependences of 
(i) the ionization rate W, which increases at shorter wavelengths and (ii) the peak intensity, 
which decreases in the UV. This decrease is due to the fact that, according to Equation 
(EJ) and considering the typical dispersion curves in gases, the higher-order indices increase 
faster, in absolute values, when the frequency increases. As a consequence, the clamping 
intensity canceling the Kerr effect iAn Kerr =o is lower, as shown on Figure |3j 

The spectral dependence of the peak electron density by more than one order of magni- 
tude, also evidenced on Figure Ufa), appears in contradiction with numerical results obtained 



without considering the higher-order Kerr terms [31], which predict almost similar electron 
densities, close to 10 16 cm -3 , at the three investigated wavelengths (248, 800 and 1550 nm). 
Since the latter results were obtained with longer pulses (127 fs) of slightly lower power than 
in our work, a direct quantitative comparison cannot be performed. However, we expect 
that the qualitative difference between our results and those of Ref. 3l| are due to the 
inclusion of the higher-order Kerr terms in our model, and thus illustrate the influence of 
the latter. 

The very similar behavior, except for the plasma density, of two gases with comparable 
values of the non-linear refractive indices but significantly different ionization potentials, 
suggest that filamentation is driven by the hig her-order Kerr terms rather than by plasma, 



not only at 800 nm as evidenced recently 



12], but also on the whole investigated spectral 



range in both air and argon. Indeed, the on-axis £ ratio (Figure [Tfg),(h)) exceeds 1 in 
all considered conditions. Furthermore, the higher-order Kerr terms strongly dominate 
the filamentation dynamics (£ > 10) at all wavelengths above A ~400 nm in argon and 
above A ~600 nm in air, although, due to the thresholds induced by the steps in photon 
numbers required to ionize oxygen (See Equation ffTUl) ). the evolution of this behavior is not 
strictly monotonic. Due to the domination of the higher-order Kerr terms, the latter will 
govern the intensity within the filaments, so that the ionization is mostly decoupled from 
the filamentation dynamics. 

It should be noted that Miller formulae are sometimes considered to underestimate dis- 
persion 



32] . However, such correction to the dispersion curve of higher-order Kerr terms 
would result in a larger absolute value of the high-order Kerr indices in the UV. In this 
case, /An Korr =o w iH be even lower in the UV. In other words, the curves of Figure [3] will 
decrease faster on the left side of the graph. As a consequence, the higher-order Kerr indices 
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will have an even more dominant contribution to the self-guiding, as compared to that of 
plasma. Such larger dispersion would therefore reinforce qualitatively our conclusion about 
the domination of defocusing higher-order Kerr terms over the plasma defocusing. 

This new vision of filamentation provides straightforward interpretation to experimental 
observations that the spectral broadening mostly occurs at the beginning of the filaments, 
while the spectrum is little affected in the main region of self-guided propagation [33| . This 
is due to the fact that, in the self-guided propagation, the intensity is clamped close to 
-^An Kcrr =o, so that the spectral counterpart of this Kerr effect, self-phase modulation (SPM), 
is also waved out. The earlier onset of filamentation for longer wavelengths is compatible 
with the Marburger formula 34J, which predicts that the non-linear focal length describing 
self-focusing is inversely proportional to the wavelength. Furthermore, the longer filaments 
in the infrared appears to stem from the lower peak electron density (Figure HJ^a)). 

The very smooth dispersion curve of both air and argon above 500 nm 2l| results in a very 
slow variation of this clamping intensity /An Ke rr=o (Figure E]) and explains the quasi-constant 
clamping intensity observed over this spectral range. 

The availability, for the first time, of simulations at constant reduced power P/P cr over 
a broad spectral range allows to discuss the choice of wavelength to optimize the filaments 
properties for specific purposes. As described above, the spectral dependence of the electron 
generation contrasts strongly with the common expectation that, due to a more efficient 
ionization, ultraviolet wavelengths should ionize the propagation medium more efficiently 

However, the peak electron density is not directly relevant for typical atmospheric appli- 
cations. For example, laser-assisted water condensation 35[, can be expected to require the 
largest possible value of the total generated charge. This total charge is obtained by inte- 
grating the plasma density over the filament length and cross-section (Figure H^b)). Here, 
the longer filaments in the infrared as well as their larger diameter result in an unanticipated 
larger total charge at longer wavelengths. 

Applications such as lightning control 36[ require the longest possible filaments with 
the higher conductivity, i.e. a high longitudinal (transversely integrated) electron density. 
Long wavelengths simultaneously optimize both of these properties (Figure HIc,d)). In that 
regard, UV lasers providing short filaments with a small diameter and decreasing transverse- 
integrated electron densities would expectedly be less efficient than the commonly used Ti:Sa 
lasers around 800 nm, while wavelengths further in the IR should be even more efficient. 
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FIG. 4: Wavelength dependence of ionization in laser filamentation. (a) Peak electron density; 
(b) total number of electrons; (c,d) Transverse-integrated longitudinal electron density in air and 
argon, representative of the electric conductivity of the filament. Note the difference with Figures 
DJc,d), which is due to the larger filament diameter in the infrared. 

Infrared is also very attractive because it would alllow to work in the "telecom" spectral 
region (1.55 /xm), where optical components are available and eye-safety standards are much 
more favorable for atmospheric applications than in the rest of the NIR and visible spectral 
region. 

In conclusion, applying generalized Miller formulas to estimate arbitrary-order non- 
linearities from the infrared to the ultraviolet spectrum, we have shown that higher-order 
nonlinear indices (up to ns in air and nio in argon) have a dominant contribution to both the 
focusing and defocusing terms implied in the self-guiding of 30 fs laser pulses in air, at all 
wavelengths between 300 and 1600 nm. As a consequence, the plasma generation is almost 
decoupled from the self-guiding of filaments over the whole visible and infrared spectral 
range. Instead, filaments can be considered as a self-guiding regime mostly governed in air 
and argon by the balance between the alternate signs of the nonlinear indices, respectively 
resulting in Kerr self-focusing and self-defocusing. 
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Moreover, a systematic investigation at constant reduced power P/P cr as a function of 
wavelength provides for the first time hints to choose optimal wavelengths for generating 
filaments optimized for a specific applications. In particular, we have shown that, contrary 
to previous expectations, the infrared is more favorable than the commonly used 800 nm 
wavelength if long ionized filaments or high total amount of charges are required. This 
raises the hope to further improve the spectacular results jl], 0, 36] recently obtained with 
Titanium: Sapphire lasers. Moreover, the higher eye-safety threshold in the telecom window 
at 1.55 fim is favorable for the practical development of the evisioned applications in free 
space. 
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